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We present a simple one-dimensional quantum-mechanical model for a particle attached to a 
surface. We solve the Schrodinger equation in terms of Weber functions and discuss the behavior 
of the eigenvalues and eigenfunctions. We derive the virial theorem and other exact relationships 
as well as the asymptotic behaviour of the eigenvalues. We calculate the zero-point energy for 
model parameters corresponding to H adsorbed on Pd(100) and also outline the application of the 
Rayleigh-Ritz variational method. 

PACS numbers: 03.65.Ge, 34.35. +a 

I. INTRODUCTION 



In an introductory course on quantum theory one commonly discusses some of the simplest models, such as, for 
example: free particle, particle in a box, harmonic oscillator, tunnelling through a square barrier, etc. with the 
purpose of making the students more familiar with the principles or postulates of quantum theory. Some time ago 
Gibbsi introduced the quantum bouncer as a model for the pedagogical discussion of some of the relevant features of 
^ ' quantum theory. He derived the solutions to the Schrodinger equation in terms of Airy functions and obtained the 
energy spectrum for two different model settings. 
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, and later by Mei and Lee^. This model is as simple as the quantum bouncer and both can therefore be discussed in the 
same course. An interesting feature of this model is that it may in principle be useful to simulate a particle attached 

■ to a wall like, for example, an atom adsorbed on a solid surface^. Although it is an oversimplified one-dimensional 
< ^ ' model of the actual physical phenomenon we deem it worthwhile to discuss some of its properties in this paper. 

In Sec. [TT] we introduce the model, write the Schrodinger equation in dimensionless form, and discuss some of 
the properties of its solutions. In Sec. IIIII we obtain the eigenfunctions and eigenvalues explicitly in terms of the 
well known Weber functions^ and show the behavior of the eigenvalues and excitation energies with respect to the 
distance between the particle and the wall. We also calculate the zero-point energy for values of the model parameters 
corresponding to de adsorption of H on Pd(100)^ and outline the application of the Rayleigh-Ritz variational method 
to the Schrodinger equation. Finally, is Sec. IIVI we give further reasons why this model may be useful in a course on 
quantum theory. 



II. SIMPLE MODEL 



We consider a simple one-dimensional model for a particle of mass m attached to a surface located at x = that 
separates free space [x > 0) from the bulk of the material (x < 0). Therefore, we assume that the potential exhibits 
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an attractive tail for x > that keeps the particle in the neighborhood of the surface and a repulsive one for x < 
that prevents the particle from penetrating too deep into the material. Since we want to keep the model as simple as 
possible we choose 

oo if x < 

V(x) = { 1 
1 §(z-d) 2 if x > 

where k 7 d> 0. We note that the particle oscillates about x = d but the motion is not harmonic because of the effect 
of the hard wall. 

The Schrodingcr equation reads 

- ^"{x) + \{x-dfi>(x)=Ei>{x) 

ip{0) = 0, lim ip(x) = (2) 

x— >oo 

where the boundary condition at x — comes from the fact that in this simple model the particle cannot penetrate 
into the material and, consequently, ip( x ) = for all x < 0. Since it is more convenient to work with a dimensionlcss 
equation we define the length unit L = I J and the dimensionless coordinate q = (x — d)/L. Thus, the 
dimensionless Schrodinger equation reads 

- \<p"{<i) + ^VO) = «p(q) 

ip(-qo) = 0, lim tp(q) = (3) 



where 
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qa = L = d Kw 

<p(q) - i>(Lq + d) (4) 

It is worth noting that qo increases with d, m, and k. It is precisely Eq. ([3]) that was discussed by Dearth and Mei 
and Lee^. 

When (/o — > oo we have the well-known harmonic oscillator with eigenvalues 

lim e n = n + -, n = 0, 1, . . . (5) 

On the other hand, when go = we have the harmonic oscillator in the half line and 
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lim e„ = 2n + -, n = 0, 1, . . . (6) 

Note that these are merely the harmonic-oscillator eigenvalues with odd quantum number (the corresponding eigen- 
functions have a node at origin). It follows from Eq. (|A.6j) in the Appendix that 

de , de 

— < 0, lim — = (7) 

dqo 90^00 dqo 

from which we conclude that the energy eigenvalues decrease monotonously between the following limits: 

3 1 

2n+->e„(g )>n + -,n = 0,l,... (8) 
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when < go < oo. It should be kept in mind that the number of zeros of y> n (q) remains unchanged as q$ goes from 
to oo. 

Equation ()A.6|) is useful for obtaining the asymptotic behavior of the energy as go oo. To this end we simply 
integrate it from — oo to —go: 

1 f~ qo w'(b) 2 

e -<», = e -oo) + - / r J\L, db (9) 

For the ground state we expect that <p(q) — > Ne~ q I 2 as g — > oo. The normalization factor is approximately given by 

N- 2 (qo) = r CO 2 dq = ^[l + e r% )] (10) 

where erf(z) is the error function. Since erf(z) < crf(oo) = 1 we write erf(g ) = 1 — £ and expand N 2 in a Taylor 
series about £ = 0: 

Art ^2_2 + £ 3 - erf ( g0 ) 

N{qo) W 20F = (U) 

Finally, Eq. © yields 

..<-») 

In order to obtain this result we substituted erf (go) ~ 1 after the integration and neglected a term proportional to 
e~ 2q ° because it is much smaller than the exponential one retained in Eq. (|T2|). This result agrees with the one derived 
by Mei and Lee^ by means of perturbation theory (note that their parameter R is V^Zo)- 
Proceeding in the same way for the first excited state we obtain 

*<-*>" I + ag T7T £j < 13) 

that is slightly different from the result of Mei and Leei However, they are equivalent for most purposes because the 
difference between them is smaller than their absolute errors. 
Eq. (|A.10p gives us the virial theorem^ for this model 

*> + <!■>—£ M 

where D = d/dq. The right-hand side of this equation is the virial of the force exerted by the surface. Eq. (|A.13[) 
provides us with another interesting relation 

<* " "li (15) 

that clearly reveals the asymmetry of the interaction between the particle and the surface. In both cases we recover 
the well-known results for the harmonic oscillator \D 2 ) + (g 2 ) =0 and (g) = when go — > oo. 



III. RESULTS 



If we define the new independent variable z = \J~2q and write the energy as e = m + 1/2 then we realize that 
D m (z) = (p(z/\/2) is a solution to the Weber equation^ 

D'; n (z) + (m + 1/2 - z 2 /A) D m (z) = (16) 
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The general solution is£ 

D m (z) = 2™/ 2 0Fe-* 2 / 4 

where the confluent hypergeometric function F (a \c\ z) is a solution to 

zF"(z) + (c - z)F(z) - aF(z) = 

and can be expanded in a Taylor series about z = as 

F (o cz ) = l + -z+ ^ ^ > z * + 
c 2!c(c + f ) 



\2z ( f — m 
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The boundary condition tp{— go) 
For each value of 90 we solve 

f 



enables us to calculate the eigenvalues from the roots of D Ti 



r( 



1 — m 
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) V 2 



2g 



r(-f) 



F 







(17) 

(18) 
(19) 

V2g ) = 0. 
(20) 



for m and then calculate the dimensionless energy e = m + 1/2. This approach has already been discussed by Dean^ 
and Mei and Lee^. 

Fig. Q] shows e n (qo) for n = 0,1,2,3 and a wide range of values of go- We note that the dimensionless energy 
decreases monotonously as predicted by Eq. ([7]) between the limits indicated in Eq. (|8|). 

The gap between two consecutive energy levels of the harmonic oscillator is Ae„ = e n +i — e n = 1 (A.E n = Hlu). 
Fig. [2] shows Ae n (qo) for n = 0, 1, 2 where we note that the energy gap increases with n revealing that the presence 
of the wall results in an anharmonic oscillation. It becomes more harmonic as qo increases (by increasing either d, m 
or k). 

In order to have a clearer physical idea of the kind of predictions of this simple model we may choose the parameters 
for the adsorption of hydrogen on Pd(100). Gladys et al^ estimated d pa 0.4 A and k « 15 N m" 1 for hydrogen that 
lead to go ~ 1-55. The zero-point energy for such model parameters is approximately 0.57 hoj instead of the value 
Tiuj/2 chosen by those authors. Although they fitted the potential energy of the vertical displacement of the H atom 
from the Pd surface to a cubic polynomial they simply chose the zero-point energy of the harmonic oscillator. Present 
model predicts that the zero-point energy is slightly greater than the harmonic-oscillator one because of the repulsive 
effect of the surface. We agree that a hard wall may not be the most adequate representation of the short-range 
interaction between the H atom and the Pd surface, but we think that the model is a reasonably simple first approach 
to the physical phenomenon. 

For deuterium we have d ~ 0.45 A which, together with the same force constant and about twice the mass, yields 
go ~ 2 and eo « 0.52 fou that is closer to the harmonic-oscillator zero-point energy^. We note the effect of the 
distance to the surface and the mass of the particle on the vibrational energies. 

This model is also useful for discussing the Rayleigh-Ritz variational method^. If, for example, we choose the trial 
function 

JV-l 

«(?) = E c ^(«) ( 21 ) 
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as a linear combination of the non orthogonal basis set 

fi(<l) = (<l + <lvWe~ q2/2 , .7=0,1,... (22) 



and minimize the approximate energy 



■w 



f_ qo u{q)Hu{q) dq 

I- qo u (i) 2 d <i 

2dq 2 ' 2' 



1 d 2 1 

H=-14^ + W (23) 



then we arrive at the secular equations 

N-l 



{Hij - wSij) Cj = 0, i = 0, 1, . . . , N - 1 

3=0 



H l3 = Mq)Hf J (q)dq, S tj = Mq)fj(q)dq (24) 

J—qo J-qo 

There are nontrivial solutions for the coefficients Cj only if the secular determinant vanishes 

\Hii-wS ij \?7j = (25) 

Its roots , n = 0, 1, . . . , N — 1 are upper bounds to the actual eigenvalues and satisfy 

«W > > e n (26) 

The particle attached to a wall is therefore a suitable example for illustrating how the roots of the secular determinant 
approach the eigenvalues (calculated accurately from the Weber functions) from above as N increases. We do not 
show results here and simply mention that the calculation is greatly facilitated by the fact that one calculates the 
integrals and SV,- analytically. 

IV. FURTHER COMMENTS 

The model discussed here is suitable for a course on quantum theory because it does not require much more 
mathematical background than it is necessary for the discussion of the well known harmonic oscillator or the quantum 
bouncer—. It is useful for introducing a numerical calculation of the eigenvalues that the student does not find in the 
treatment of the harmonic oscillator. One can approach the problem by means of either the Weber functions or the 
Rayleigh-Ritz variational method. The student will also learn that it is necessary to modify the form of the well 
known virial theorem and other mathematical expressions in order to take into account the effect of the wall. We 
believe that present derivation of the analytical results in Sec. [IT] and i n the Appendix is simpler than those available 
in the scientific literature^. 

In addition to it, the model enables us to simulate the adsorption of an atom on a surface and discuss anharmonic 
vibrations in quantum theory. In the study of molecular vibrations one introduces nonlinear oscillations by means 
of cubic, quartic, and other terms of greater degree in the potential-energy function. In this case it arises from the 
boundary condition forced by the hard wall. 



Appendix: Some useful mathematical relations 

In this appendix we develop some useful analytical results for the eigenfunctions and eigenvalues of the constrained 
oscillator. Although similar expressions have already been shown elsewhere^ we derive them here in a form that is 
more suitable for our needs. 

Consider the dimensionlcss Schrodinger equation 

- \f"i x ) + V(x)<p(x) = etp(x) 

<p(b) = 0, lim <p(x) = (A.l) 

Both the eigenvalue e and the cigenfunction ip(x) depend on the chosen value of b. If we differentiate this equation 
with respect to b and call x( x ) = dip(x)/db we have 

- \x"(x) + V(x) X (x) = e X {x) + ^ip(x) (A.2) 

If we multiply this equation by ip(x) and integrate the result between b and oo we easily obtain 

- J b V {x) 2 dx = -- X (b) v '(b) (A3) 

because the integration by parts of the first term yields 

x"(x)tp(x)dx = X (b)v'(b) + 2 [V(x)-e) X (xMx)dx (A.4) 

b Jb 

If we now differentiate the boundary condition (p(b) = and take into account that f(x) depends on b we obtain 

X (b) + <p'(b)=0 (A.5) 



so that equation (|A.3[) becomes 



We define the operators 



| / V{x?dx=-v'{bf (A.6) 



H = -\b 2 + V{x) (A.7) 

and v = xD, where D = d/dx. The commutator between them reads 

[H, v]=Hv- vH = -D 2 - xV (A.8) 

If tp(x) is an eigenfunction of H with eigenvalue e then straightforward integration by parts leads to 

ip[H,v](pdx = --bip'{bf (A.9) 
6 2 



which, by virtue of Eq. (IA.6I) , becomes the virial theorem 



+ (XV) = b^ (A.10) 



where 



- - (A.11) 
Jb V dx 
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Analogously, from the commutator [H,D] = — V wc obtain 

i 

y{x) 2 V'{x)dx = -<p'{b) 2 (A.12) 



(V) = | (A.13) 

We can test these equations quite easily by means of the well-known solutions to the free harmonic oscillator 

%l> n = N n H n (x)e- x2 / 2 (A.14) 

where H n (x) is a Hermite polynomial and N n the corresponding normalization factor—. If we choose b to be one of 
the zeroes x n j, j = 1,2, ... ,n of H n {x), n = 1,2,... then wc verify that i/>„(x) satisfies equations (|A.9[) and (|A.12[) 
for V(x) — x 2 /2. If, for example, b = for n = 1 and b = ±1/V2 for n = 2 then e = 3/2 is the energy of the ground 
state of the harmonic oscillator with the boundary condition at b — and e = 5/2 is the energy of the ground state 
with b = \j\pl and of the first excited state with b = —l/y/2. 
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FIG. 1: (Color online) First four energy levels 




FIG. 2: (Color online) The first three excitation 6I16rgl6S €n -j-l 671 (n — 0, 1, 2 from bottom to top) 



